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Molecular diffusion between walls with adsorption and desorption 
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The time dependency of the diffusion coefficient of particles in porous media is an efficient probe of their 
geometry. The analysis of this quantity, measured e.g. by nuclear magnetic resonance (PGSE-NMR) , can 
provide rich information pertaining to porosity, pore size distribution, permeability and surface-to-volume 
ratio of porous materials. Nevertheless, in numerous if not all practical situations, transport is confined by 
walls where adsorption and desorption processes may occur. In this article, we derive explicitly the expression 
of the time-dependent diffusion coefBcient between two confining walls in the presence of adsorption and 
desorption. We show that they strongly modify the time-dependency of the diffusion coefBcient, even in this 
simple geometry. We finally propose several applications, from sorption rates measurements to the use as a 
reference for numerical implementations for more complex geometries. 



I. INTRODUCTION 

The time dependency of the diffusion coefficient of 
particles in porous media is an efficient probe of their 
geometry^^-. Extracting information on the microscopic 
structure of the porous medium from diffusion measure- 
ments, has seen applications from heterogeneous geo- 
logical materials, e.g. clays^ to biological tissues, e.g. 
brainS. In this respect, pulsed gradient spin echo nu- 
clear magnetic resonance measurements (PGSE-NMR), 
unraveled by explicit derivations of the time-dependent 
diffusion coefficient of tracers in simple, model media, "^'^'^ 
have shown to reveal rich information regarding to poros- 
ity, pore size distribution, permeability and surface-to- 
volume ratio.— 

Nevertheless, in numerous if not all practical situa- 
tions, transport is confined by walls where adsorption 
and desorption processes may occur. This is especially 
true for porous media for which the surface-to- volume ra- 
tio is large. Such processes dramatically alter the overall 
dynamics of the diffusing species, even in the most simple 
geometries^iiSrJi. 

Here, we derive explicitly the expression of the time- 
dependent diffusion coefficient between two confining 
walls in the presence of general adsorption and desorp- 
tion. 



II. RESULTS AND DISCUSSION 

Following a stochastic approach, we consider a Brown- 
ian particle moving along the x axis with a bulk diffusion 



coefficient Db, confined between two walls at positions 
X = and x = L. The evolution equations and the 
boundary conditions are 



dtc{x,t) = Db-g^ for X e ]0,L[, 

hit) = -kdTo + kaC (0, t) = Db dxc\^^o , 
rL(t) = -kdTL + kaC (L, t) = -Db d^cl^ 



L ' 



(1) 

(2) 
(3) 



where c{x,t) [length^'^] is the particle concentration 
at position x and time t, Tq and [length^^] are 
the surface concentrations in adsorbed particles at both 
walls. The desorption and adsorption rates are noted kd 
[time~^] and ka [length-time^^]. The brownian particles 
can be in three states. In state ab, they're free to move 
along the x axis. In states oq or a^, they're adsorbed 
at one wall at a; = or L, where they remain immo- 
bile. The equilibrium probability Peq associated to these 
states follows from the normalization condition: 



^ i cq 1 

=ao,aL 



Pcq (a) + / -Pcq (afo) dx = 1, 



(4) 



with the result Pcq (do) = ^cq (ol) = ka/{kdL + 2ka) and 
-Peq (a&) = kd/{kdL + 2ka). 

Assuming that the particles are initially distributed 
according to Peq, the mean squared displacement (MSD), 

M{t) = (^{x (t) — xq)"^^, of the particles can be expressed 
as 



J 



M(i) 



dx 



dxo [x - xo f ^ Peq (a) G [x, /3, t\xo, a) . 

a, (I 



(5) 
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where G {x, /3, t\xQ, a) is the particle propagator, i.e. the 



2 



probability for a particle initially at Xq in state a to be 
at X in state /3 at time t. Several propagators have to 
be calculated, which correspond to all combinations of 
starting and ending positions and states. Their explicit 
determination is conveniently performed by first Laplace 



transforming Eq. HHSl In the Laplace domain, bound- 
ary conditions associated with adsorption and desorp- 
tion simplify to radiative boundary conditions (see, for 
example, KeiJ^^). Then, after standard calculations, 
the Laplace transform of these propagators are found to 
be given by: 



G{x < xo,ab,s\xo,ab) = 
G{x > Xq, ab, s\xo,ab) = 
G {x, ao, s\xo,ab) = 


Kq{j cosh((7x) + q ainh{qx)) (7 cosh (q {xq — L)) - 
Kq (7 cosh (qxa) + gsinh (qxo)) (7 cosh(g(a; — L)) 
—q {jka — s) (7 cosh{q{L — x j) + q sinh(g(L — x 


g sinh (g (.To - 
— q sai\i{q{x — 

))), 


m, 


(6) 
(7) 
(8) 


G {x,aL,s\xo,ab) = 


—q{"fka — s) (7C0sh(ga;) + g'sinh(ga;)). 






(9) 


G {x,ab,s\xo,ao) = 


-^q {ika - q^Db) (7cosh(g(xo - L)) ~ qsu\\\{q 


(^o-i))), 




(10) 


G {x, ao, s\xo, ao) = 


-r-TTl (7^a - s) {{-fka + q^Db) cosh{Lq) + q {ka 


^- 7Db) sinh(Lg)) , 


(11) 


G {x,aL,s\xo,ao) = 


, , {jka s) {-fqka q^Db) , 






(12) 


G {x,ab,s\xQ,aL) = 


{l^a ~ q^Db) (7 cosh {qxo) + q sinh {qxo)) , 

l^d 






(13) 


G {x,aQ,s\xo,aL) = 
G {x,aL,s\xo,aL) = 


G {x,aL,t\xfi,ao) , 
G(x, ao,i|xo,ao) , 






(14) 
(15) 



where s is the conjugate of t, q = y^I/D^, 7 ^ ^^-^ and k ^ {s (7^ + q^) sinh(Lg) + 2sjq cosh(Lg)) ^ . 

ka Inserting these propagators and their respective Poq in 

Eq. O one gets for the Laplace transform of the mean 
squared displacement: 



M{s) 



1 


-2kdL 4fcd(fcd + s)sinh 


m 




X 

2fca + k(iL 


- ■''1^ sq3 ((fcrf + s)cosh 


+ kaq sinh(^)) 



(16) 



where s is the conjugate variable of time t in Laplace 
space, and q — \J-^- The diffusion coefficient is then 
easily derived from the MSD as 

Dit) ^ (^^Mis)) , (17) 

where is the inverse Laplace transform. The result- 
ing reduced time-dependent diffusion coefficient is plot- 
ted in figure [T] as a function of the reduced time Dbt/L^, 
for various adsorption and desorption rates^. At t = 0, 
D{0)/Db is simply the fraction f = {1 + 2ka/kdL)-^ 
of mobile (non-adsorbed) particles. Figure [T] explores a 
large range of /, from weak to strong sorption. For all 
sorption parameters ka and fc^, diffusivity decreases with 
time, as increases the probability for each Brownian par- 
ticle that explores the media to experience confinement 



by a wall and immobilization by adsorption. This de- 
crease is strongly modified by the sorption processes, as 
the more probable it is for a particle to get adsorbed, the 
more time it takes to explore the media. At long times, 
because our media is not permeable, i.e. because the con- 
finement is total in our geometry, the effective diffusion 
coefficient tends toward zero. 



III. SHORT AND LONG TIMES APPROXIMATIONS 

-'^^'^ FT^"*" 2FT ^ 1' ^^-"^ sufficiently long residence 
time on walls and fraction of adsorbed species, the series 
expansion of the hyperbolic functions in Eq. [TC] allows, 
together with the residue theorem, to explicitely inverse 



3 



1.0 



0.4 



0.2 - 



0.0 





1 


/ = 10-^ : 










\\ 




\\ 










0.5 


/=^^^^__ 





0.0 



0.2 



0.4 

Db-t/L 



0.6 

2 



0.8 



1.0 



FIG. 1. Time-dependent diffusion coefBcient in reduced units 
as a function of the reduced time. Several ratios / = (1 + 
2ka/kdL)~^ are reported, which decrease with the strength 
of the sorption phenomena: (black line) / = 1, i.e. without 
adsorption; (red dashed line) / — 0.5; (green dot-dashed line) 
/ = 10"^ and (blue dotted line) / = 10"^ 
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FIG. 2. Reduced diffusion coefficient as a function of reduced 
time for / = 10~^; (black) from the full expression of Eg. 1171 
(red dashed line) from the short time approximation of Eg. 1191 
up to order 10 and (blue dotted line) from the long time ap- 
proximation of Eq. 1181 



Eq. [T7]in the long time-limit {i.e. for s — > 0), leading to 

kd 



kaL 



1 



kgL 
215. 



t 



(18) 



At short times, we analytically inverse the power series 
expansion for C^^ [lims_5.oo D (s)]. The first terms read: 



D. 



short 



it) 



= /x 




(19) 



See Supplementary Material Document Nojii for terms 
up to order 10 of the expansion. Figure [5] compares 
the short and long time approximations of the time- 
dependent diffusion coefRcient to the full numerical inver- 
sion of the Laplace transform of Eq.[T7]for / = 10~^. To- 
gether, the short and long time approximations provide 
a simple uniform approximation over the whole temporal 
domain. 



IV. CONCLUSION 

In summary, we addressed explicitely the effect of 
adsorption and desorption phenomena on the time- 
dependent diffusion coefficient of Brownian particles con- 
fined between walls. The full expression of D{t) is given, 
together with convenient and accurate approximations. 
It is shown that adsorption and desorption processes 
strongly modify -D(t), so that their consideration is im- 
portant for proper determination of micro-geometric in- 
formation from D(t) measurements. Eg. [TTlmav also be 
used to extract sorption rates ka and kd from experimen- 
tal measurements of D{t), e.g. by PGSE-NMR. Finally, 
in more complex geometries, resort to numerical simu- 
lations such as Lattice Boltzman n^^i^^ is necessary. The 



present work provides exact reference results required for 
the validation of such numerical schemes in the presence 
of adsorption and desorption. 
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I. SHORT-TIME EXPANSION OF THE TIME-DEPENDENT DIFFUSION 
COEFFICIENT 



At short times t, the power series expansion for the time-dependent diffusion coefficient 
up to order 10 reads: 



■Dshort(*) _^ 4.^Vt 2kat 

f ■ Di LxAF L 

kg {Djka - fcg) 
16fcg (2D,fcd - fcg) 

ka {kt - SDtkakl + Dlkfj 
32fe2 (Di,fcd - kl) {sD^ka - kl) 
kg (-fcg + bDbkgki - 6Dgfcgfcg + ggfcg) t4 
^ 64fcg (2i3fcfcd - fcg) (fc^ - 4ij;,fcdfcg + 2Dgfcg) t^/^ 

fcg {Pbkd - fcg) (-fcg + GDbkgkt - Qggfcgfcg + ggfcg) 

60LD4 

128fcg (fc^ - SD^kdkl + ggfcg) (fcf _ 5Dikdkl + SJgfcg) f"/" 
103951, VttD^/^ 

fcg (-C + ggfefcdfcg - 28Dgfcgfcg + 35ggfc3fc^ - ISPgfc^fcg + ggfeg) f6 

360LD5 

^ 256fcg (gfefcd - fcg) (2Dj,fcd - fcg) (sDfcfcrf - fcg) (fc^ - 4Di,fcrffcg + ggfcg) t^^/^ 

135135L0FdJ^''^ 

fcg (fcl" - llJfefcdfc^" + 45J'gfcgfcg - 84g3fcgfcg + 70J'4fc4fc4 _ 21g5fc5fc2 ^ j,6fc6^ ^7 

2520LD6 

512fcg (-fcg + bPbkqkt - 6Dgfcgfcg + Dgfcg) (-fcg + 7Dbkakt - 14Dgfcgfcg + TDgfcg) t^^/^ 

2027025Lv^d"/^ 

fcg (gfefcd - fcg) {kt - 3D6fcdfcg + ggfcg) (fcg - 9gbfcdfcg + 26Dlklkt - 'ilDlklkl + g^fc^) 

20imLDl 

1024fcg (2J),fcd - fcg) {kt - iD.kgkl + 2£>gfc2) (fcg - SDtfcdfcg + 20J'gfcgfcf - 16ggfcgfcg + 2D4fc4^ ^17/2 

34459425LV^I?^^/^ 

fcg (fcge - ISJfcfcrffcM + 91ggfcgfcg2 - 286D3fc3fcg° + 495g4fc4fc8 _ 462D5fc5fcg + 210Dgfcgfc4 - 36g^fc^fcg + Dgfcg) 

181440LDf 

2048fcg (P.fcd - fcg) {3D,ka - fcg) (-fcg + GD.kgkt ~ 9ggfcgfcg + Pgfcg) (-fcg + 6J'fcfcrffc4 - QZPgfcgfcg + 3J'3fc3^ ^19/2 

654729O75L0rZ)"/^ 

+0(tl°), 
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where Db is the bulk diffusion coefficient, ka and kd are the adsorption and desorption rates, 
L is the distance between the confining walls, and / = (1 + 2ka/kdL)'~^. 
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